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The construction of the semiclassical trace formula for the resonances with the transverse electric 
(TE) polarization for two-dimensional dielectric cavities is discussed. Special attention is given to 
the derivation of the two first terms of Weyl's series for the average number of such resonances. The 
obtained formulas agree well with numerical calculations for dielectric cavities of different shapes. 
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I. INTRODUCTION 

Open dielectric cavities have attracted a large interest 
in recent years due to their numerous and potentially 
important applications [l|, [1] . From a theoretical point of 
view, the crucial difference between dielectric cavities and 
much more investigated case of closed quantum billiards 
that in the latter the spectrum is discrete but in 
the former it is continuous. Indeed, the main subject of 
investigations in open systems is not the true spectrum 
but the spectrum of resonances defined as poles of the 
scattering S'-matrix (see e.g. d, 01)- 

The wavelength of electromagnetic field is usually 
much smaller than any characteristic cavity size (except 
its height) and semiclassical techniques are useful and 
adequate for a theoretical approach to such objects. It is 
well known that the trace formulas are a very powerful 
tool in the semiclassical description of closed systems, see 
e.g. Therefore, the generalization of trace formu- 

las to different open systems, in particular to dielectric 
cavities, is of importance. 

The trace formula for resonances with transverse mag- 
netic (TM) polarization in two-dimensional (2d) dielec- 
tric cavities has been developed in and shown to agree 
well with the experiments and numerical calculations 

This paper is devoted to the construction of the trace 
formula for 2d dielectric cavities but for transverse elec- 
tric (TE) polarization. Due to different boundary con- 
ditions the case of TE modes differs in many aspects 
from TM modes. In particular, a special treatment is re- 
quired for the resonances related to Brewster's angle [ll| 
at which the Fresnel reflection coefficient vanishes. 

Our main result is the asymptotic formula in the semi- 
classical (aka short wave length) regime for the average 
number of TE resonances for a 2d dielectric cavity with 
refraction index n, area A and perimeter C 
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Here iVTE(^) is the mean number of resonances (defined 
below) whose real part is less than fc, the coefficient tte 



is given by the expression 
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and i?TE is the Fresnel reflection coefficient for the scat- 
tering on a straight dielectric interface at imaginary mo- 
mentum 



RTEit) = 
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The plan of the paper is the following. In Sec. |TI] the 
main equations describing the TE modes are reminded. 
In Sec. lIIII the circular cavity is briefly reviewed: an exact 
quantization condition is derived, which allows a direct 
semiclassical treatment. In Sec. IIVI the first two Weyl 
terms for the resonance counting function are derived. It 
is important to notice that, for TE modes, one can have 
total transmission of a ray when the incidence angle is 
equal to Brewster's angle. This leads to a special set of 
resonances, which are counted separately in Sec.|Vl Sec- 
tion IVII is devoted to a brief derivation of the oscillating 
part of the resonance density. In Sec. IVIII our obtained 
formulae are shown to agree well with numerical com- 
putation for cavities of different shapes. In Appendix |X] 
another method of deriving the Weyl series for TE po- 
larization based on Krein's spectral shift formula is pre- 
sented. 



II. GENERALITIES 

To describe a dielectric cavity correctly one should 
solve the 3-dimensional Maxwell equations. In many ap- 
plications the transverse height of a cavity, say along the 
z axis, is much smaller than any other cavity dimensions. 
In such situation the 3-dimensional problem in a rea- 
sonable approximation can be reduced to two 2d scalar 
problems (for each polarization of the field) following the 
so-called effective index approximation, see e.g. ITsj 
for more details. 
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In the simplest setting, when one ignores the depen- 
dence of the effective index on frequency, such 2d ap- 
proximation consists in using the Maxwell equations for 
an infinite cylinder. It is well known [ll'] that in this ge- 
ometry the Maxwell equations are reduced to two scalar 
Helmholtz equations inside and outside the cavity 



{A + n^k^)-i'{x) = 0, 
{A + k^)'^{x) = 0, 



xeV, 
X 4V , 



(4) 



where n is the refractive index of the cavity, T) indicates 
the interior of the dielectric cavity, and — for the 
TM polarization and ^ = for the TE polarization. 

Helmholtz equations ^ have to be completed by the 
boundary conditions. The field, ^'(a?), is continuous 
across the cavity boundary and its normal derivatives 
along both sides of the boundary are related for two po- 
larizations as below [ll| 
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Open cavities have no true discrete spectrum. Instead, 
we are interested in the discrete resonance spectrum, 
which is defined as the (complex) poles of the ^-matrix 
for the scattering on a cavity (see e.g. Q)- It is well 
known that the positions of the resonances can be deter- 
mined directly by the solution of the problem (|4]) and ([5]) 
by imposing the outgoing boundary conditions at infinity 
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The set ([l])-® admit complex eigen- values k with Im k < 
0, which are the resonances of the dielectric cavity and 
are the main object of this paper. Our goal is to count 
such resonances for the TE polarization in the semiclas- 
sical regime. This will provide us with the analogue of 
Weyl's law derived for closed systems, see e.g. [T^ . 



III. CIRCULAR CAVITY 

The circular dielectric cavity is the only finite 2d cavity, 
which permits an analytical solution. Let R be the radius 
of such a cavity. Writing 'i'{x) — AJ,n{nkr)e™'^ inside 
the cavity and ^(a? ) — BHm^ (fcr)e""'^ outside the cavity, 
it is plain to check, that in order to fulfill the boundary 
conditions, it is necessary that k is determined from the 
equation Smix) = with x — kR and 



(x) ^x\-f„^{nx)Hil\x) - J™(nx)il,«'( 



(^) (7) 



where Jm{x) (resp. Hm\x)) denotes the Bessel function 
(resp. the Hankel function of the first kind). Here and 
below the prime indicates the derivative with respect to 



the argument. Factor x in ([7]) is introduced for further 
convenience. 

Using Jm(x) ~ {Hrn\x) + Hm^)/2 the equation 
Smix) = can be rewritten in the form 



Rr,r{x)Ern{x) = 1 



where 



and 
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In the semiclassical limit, a; — >■ oo, the asymptotic for- 
mula for the Hankel function [l^ {0 < m < x) gives 



H^r^Hx) 

where 



2)1/4 



5'*'"(") [l + 0(a;-i)] (11) 



'I'm (a:) — \/ x^ ~ m? 
In this way one obtains 

Em. - 
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and 



-Km (a;) 



i?TE (^) 



(12) 



(13) 



(14) 



where _Rte is the standard TE Fresnel coefficient for the 
scattering on an infinite dielectric interface 
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The above formulas mean that in the semiclassical limit 
Eq. ([U takes the form 



-Rte ( — 
V X 
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or 
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with integer p = 0, 1, 2, . . .. 

In fact, this equation is valid in the semiclassical limit 
for closed and open circular cavities with other bound- 
ary conditions as well. The only difference is that, in- 
stead of the Fresnel reflection coefficient i?TE, it is nec- 
essary to use the reflection coefficient for the problem 
under consideration. For example, for closed billiards, 
n — 1 and for Neumann (resp. Dirichlet) boundary 
conditions Rm{x) in (ITi|) equals to 1 (resp. —1). For 
open dielectric circular cavity with the TM polarization 
Rm{x) Rtm i'm/x), where Rtm is the usual Fresnel 
refiection coefficient for the TM modes [Till 



i?TM(i) 
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IV. WEYL TERMS 



Semiclassical formulas like Eq. (|T6|) are convenient to 
obtain the average number of eigenvalues and resonances 
for closed and open systems with different boundary con- 
ditions. 

Let us consider first the simplest case of a closed 
billiard with Neumann boundary conditions for which 
Rm(x) = 1. In the semiclassical regime the eigenval- 
ues for this model are determined from Eq. (jl7|) which 
reads 



np , 



(19) 



where ^m{x) is defined in and p = 0,1,... is an 
integer. Therefore, for fixed m, the number of eigenvalues 
less than x is [Njn{x)] where [ . ] stands for the integer 
part and 



n 



(20) 



1 is added as the integer p in (IT^ starts with but 
[Nm{x)] has to begin with 1. 

Summing over all m leads to the total number of eigen- 
values less than x, usually called the counting function. 
This sum is finite as the asymptotics PT|) is valid when 
\m\ < X. Finally 
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The averaged number of levels is determined from the 
equation 



(22) 



m— — x 



With a needed precision one can substitute the summa- 
tion over m by an integral and, consequently, the aver- 
aged number of eigenvalues for a circular billiard with 
Neumann boundary conditions can be approximated as 
follows 



N{x)=2 r (-{V^ 
Jo 

Using the formula 



/ (Vl - t"^ - t arccos(i))dt = ^ 
Jo o 



one gets 



N{x) 
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x + 0{l) 



— \ dm . 

4/ 
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As for the circle the area is .4 = nR^ and the perimeter is 
£ = 2itR, these results can be rewritten in the standard 
form 14 j 
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with r = 1. For Dirichlet boundary conditions similar 
arguments show that 1 /4 in (|23| is substituted by — 1/4 
and r = —1, as it should be 14]. 

For open cavities Eq. (ITB)) gives complex solutions (res- 
onances) k — ki + ik2 with negative imaginary part, 
k2 < 0. In the semiclassical limit for all investigated 
cases one has |fc2| ^ fci- Separating the imaginary and 
real parts in (fT7)) and using that 
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one gets that in the first order in fc2 the real part of the 
resonance position, ki (or xi — kiR), is determined from 
the real equation similar to (IT71) 

^„i{nxi) + S{m/xi) = irp, p integer, (28) 

where 26 (t) is the argument of the reflection coefficient 

R{t) = |i?(t)|e''*(*) . (29) 

In the same approximation the imaginary part of the 
resonance position, k2 is 



ln\R{xi)\ 



2y^ 



m? / XI 



(30) 



This semiclassical approximation is quite good even for 
not too large m as indicated in Fig. [TJ 

The above arguments demonstrate that the total num- 
ber of resonances can be calculated from the real equa- 
tion ([28|) . As in ([22]) one concludes that the mean number 
of resonances with real part xi less than x is given by the 
expression 
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Consider first the case of TM modes. The reflection co- 
efficient in this case is given by (IT51) and one has 
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FIG. 1. (Color online). The black circles are the exact positions of the resonances with m = 23 and n = 1.5 (a), n — 2 (b), 
and n = 3 (c). The blue lines indicate the approximation (|30|) . The additional levels are encircled by the red circles. The red 
stars show the approximate formula l|48[l. 



By integration by parts and contour deformation it is 
easy to check that 



2 / arctan : at = — (n — 1) 



-Rtm(^) 



1 



dt 



(35) 



where Rtm {t) is the same as but for pure imaginary 
argument 



i?TM(t) ^ i?TM(ii) = . (36) 

V'T- + <^ + V 1 + i 

Finally, these considerations lead to the expression simi- 
lar to (HH) 

^—-r / 1 \ ft k H k /-,_, 

A^TM (fc) = — + rTM -r + oik) (37) 

where 



rTM{n) = 1 + 



_ (38) 
which agrees with the result in Q obtained by a different 
method. 

Consider now TE modes. In this case the reflection 
coefficient is given by (1151) and its argument is 



arctan | — , | , 1 <t <n 



^TE(i) = < 



0, t* <t< 
0<t<t* 



1 (39) 



where t* corresponds to the zero of the TE reflection 
coefficient (Brewster's angle), -Rte(^*) = 0, 



t* 



\J-n?- + 1 



(40) 



Using these values we get 



nTV / 1 \ Tl k H k /,^\ 

^TE(fc) = — ^ +rTE-^+o{k) (41) 
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where rxE is given by the following expression 
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Similar to (|35p one can prove that 

' 2 
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, , dt = -(n- 1) 

Vn2 - i2 / 2^ ^ 



„ ^TE(i)^^2 + i2 1 + ^2 



dt 
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where, as above, i?TE(t) is the TE reflection coefficient 
()15p analytically continued to imaginary t 



~ , , X V"-^ + 12 - n2\/l + 12 

Combining all terms together we obtain that 



TTEin) = 1 + - / i?TE(t) 
Tl" J-oo 
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The first two terms are the same as for TM modes ((38)) 
but with TE reflection coefficient. The last term is the 
new one related to the change of the sign of the TE re- 
flection coefficient. 

Higher order terms in Weyl's expansions (|37| and (|4T|) 
are not yet calculated so we prefer to use a conservative 
estimate of them as o(fc) though all numerical checks sug- 
gest that for smooth boundary cavities it is 0(1). 



with integer m = 0, ±1, ±2, . . .. Therefore, the number 
of additional resonances related with Brewster's angle is 



iVadd(fc) 



Lk 



(51) 



Comparing it with Eq. ps)) we conclude that the second 
term in the Weyl expansion for the averaged number of 
resonances for TE polarization is the following 



ADDITIONAL RESONANCES 



Formula (j45|) is the correct description for the reso- 
nances whose real part of the eigen-momentum k corre- 
sponds to non-zero reflection coefficient (i.e. m/xi ^ t*). 
This is due to the fact that when the reflection coefficient 
is zero its phase is not defined. For TE modes there is 
a special branch of resonances for which semiclassically 
the real part does obey m/xi = t* . The existence of such 
additional resonances were first discussed in a different 
context in [l6j . 

The approximate positions of these resonances can be 
calculated as follows. Assume that the resonances have a 
large imaginary part. As Hm {x — it) tends to zero when 
r — > cx) one can approximate Eq. ([7]) by 

Smix) = 0, Sm{x) = —^{nx) - n—j^{x) . (46) 
From the asymptotic ([TT|) it follows that 

r(l)' 



H. 



iA 1 



2(a;^ — m^) 
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Using this expression one concludes that the solution of 
the equation Sm(im) = has the form 



1 . (n2 1)3/2 
— m — i- 



2n2 



(48) 



This approximation is better for large n when the imagi- 
nary part is large but it gives reasonable results even for 
n of the order of 1. In practice one may use (|48)) as the 
initial value for any root search algorithm (cf. Fig. [T|). 

From (^51) it follows that the ratio m/xm tends to t* 
defined in (PD|) so these resonances are not taken explic- 
itly into account in Eq. (j45p . Their number can be esti- 
mated as follows. The discussed resonances correspond 
to waves propagating along the boundary whose direc- 
tion forms an angle with the normal exactly equal to 
Brewster's angle 



smf^B = 



1 



(49) 



If the length of the boundary is £, the possible values for 
the momenta of such states in the semiclassical limit are 



kra sin 6*6 
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At 
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(52) 



where the plus sign is used when the above additional 
resonances are taken into account and the minus sign 
corresponds to the case when these resonances are ig- 
nored. 

For small values of n the additional resonances are 
mixed with other resonances and their separation seems 
artificial. For large n the additional branch of resonances 
is well separated from the main body of resonances and 
one can decide not to take them into account. In such 
a case, the minus sign has to be used in (|5^ (see below 
Section 

When the cavity remains invariant under a group of 
symmetry it is often convenient to split resonances ac- 
cording to their symmetry representations. For refiection 
symmetries it is equivalent to consider a smaller cavity 
where along parts of the boundary one has to impose 
either Dirichlet or Neumann boundary conditions. In 
this case the total boundary contribution to the average 
counting function iV(fc) is given by the general formula 



1 

47r 



{n{Lf^ - Lb) + rTEin)LQ) k . 



(53) 



Here L n and L £> are the lengths of the boundary parts 
with respectively Neumann and Dirichlet boundary con- 
ditions and Lq is the length of the true dielectric inter- 
face. It is this formula, which will be used in Section IVlII 
for dielectric cavities in the shape of a square and a sta- 
dium. 



VI. OSCILLATING PART OF THE TRACE 
FORMULA 

The quantization conditions (jT)) or ([Sj permit also to 
obtain the resonance trace formula for a circular dielec- 
tric cavity. Let kj = kij — \k2j be resonance eigen- 
momenta. Define the density of resonances as follows 
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(54) 

In general, if the zeros of a certain function Fix) 

which has no other singularities then the density of these 
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zeros ([54]) formally is given by the following expression 

(55) 



TT F{x) 

In the semiclassical limit fc — cx) it is sufficient to consider 
the semiclassical formula (ITBl) i.e. F{x) = Fm{x) and 



(56) 



A more careful discussion is performed in Appendix |^ 
In such a manner one gets 
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(57) 



Here Rm — RTEifn/x) is the Fresnel reflection coefficient 
for TE polarization (fTSj) . 

The further steps are as usual, see e.g. Q. Using the 
Poisson summation formula 



E 

m— — oo 



f{m) 
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the expression (|57|) becomes 
(i(°""'(fc) = 



/(m)dm (58) 



(59) 



where the action is 



- ^ r dmX2~^ ^Kne'^-'-(")+c.c 



SM^rim) = 2t: Mm + 2r$„i(nx) 



(60) 



When fc — )■ cx) the dominant contribution to the integral 
is due to saddle point solutions msp determined from the 
equation dSM,r{jn) / dm = 0. It is plain that 



msp ~ nxcos{6M,r) 



(61) 



with Om.t — TrM/r. This saddle point corresponds geo- 
metrically to a periodic orbit of the circle in the shape of 
regular polygon with r vertices going around the center 
M < r times. Expanding the action SM,r(jn) around the 
saddle point (|6T1) one gets 

Sm r{msp + Sm) fa nklp — —r -\ — (Jm)^ . (62) 

2 nx sm OM,r 

Here Ip — 2rRsa\9M,r is the classical length of the peri- 
odic orbit determined by M and r. 

In the end one gets the trace formula for the resonances 
of the circular dielectric cavity in the form 

s°'^\k) = 2fc— ^^=ii;:e'["^-'--'-'^/2+^/^i +C.C. 



(63) 



where Ap = ttR"^ sin^ 9M,r is the area occupied by a given 
periodic orbit family, Rp — R'x^{n cos Om.t) is the Fresnel 
reflection coefficient for the TE scattering with an angle 
equal to the reflection angle, Om.v, for the given periodic 
orbit. 

Repeating the arguments presented in Q we argue that 
in general the oscillating part of the resonance trace for- 
mula in the strong semiclassical limit has the form of the 
sum over all classical periodic orbits 



d{E) 



p 



dp{E) 



(64) 



where contribution of an individual orbit depends on the 
orbit considered 

• For an isolated primitive periodic orbit p repeated 
r times 



dp{E) 



nlr, 



TT |det(M; - 1)1 



nr irnklp 
1)11/2 -^P'- 



-ir/ip7r/2 



(65) 



where lp,Mp,^p, Rp are, respectively, the length, 
the monodromy matrix, the Maslov index, and the 
total TE Fresnel reflection coefficient for the chosen 
primitive periodic orbit. 



• For a primitive periodic orbit family 



dp{E) 



TT yj2/Kklp 



{R:)e 



inklp—irfipTc/2 



(66) 



where Ap is the area covered by one periodic orbit 
family, (i?^) is the mean value of the TE Fresnel 
reflection coefficient averaged over a periodic orbit 
family. 

The only difference with corresponding results derived in 
\§\ is that the TE reflection coefficient is used instead of 
the TM coefficient. 



VII. NUMERICAL VERIFICATION 

The numerical calculations of the resonance spectrum 
for the TE modes of the circular dielectric cavity is pre- 
sented in Fig. [2] Notice that when the cavity refraction 
index n increases the additional branch of the resonances 
(1151) separates more and more from the main part of the 
spectrum. 

In Fig. [3] we plot the difference between the func- 
tion counting the numerically computed resonances res- 
onances with a real part less than k (with radius R = 1) 
and the best fit to it of the form, see (|4T|) , 



Nfitik) = — (fci?)2 + ai(fci?) + ao 



(67) 



where ai and uq are fitting parameters. 

For n = 1.5 and n — 2 we consider all resonances 
including the additional branch. For n = 3 this branch 
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FIG. 2. (Color online). Resonance spectra for TE modes of the circular dielectric cavity with n = 1.5 (a), n = 2 (b), and n = 3 
(c). The red circles encircle the additional branch of resonances obtained by choosing initial conditions (|48|l and running a root 
searching routine to solve the equation Sm{x) = with Sm{x) given by ((7|). The range along the x axis are chosen such that 
every plot contains around 12000 resonances (counted with multiplicity). 



is quite far from the other resonances (cf. Fig. [2]c)) and 
it is natural not to include it in the counting. The fitted 
values of the parameters for these three cases are the 
following 

n = 1.5, ai = 1.246, ao = -0.66 

n = 2, ai = 1.122, ao = -0.50 . (68) 

n = 3, ai = -1.189, ao = 0.12 

The term oi has to be compared with the theoretical 
prediction which follows from ((5^ (used with plus sign 
for n = 1.5 and n = 2, and with minus sign for n = 3) 

n ^ 1.5, a*i'' = 1.247 

n = 2, af" = 1.124 . (69) 

n = 3, af = -1.190 

The agreement with our numerical calculations is very 
good. 

In Fig. the Fourier transform of the resonance den- 
sity for the circular dielectric cavity with different values 
of the refractive index is displayed. As expected from 
the trace formula, this quantity has peaks at the length 
of classical periodic orbits of the circle. Notice especially 
that the triangular orbit is not confined for n = 1.5. 
Hence the Fresnel reflection coefficient is small and in- 
duces damping, which can be clearly seen in Fig. |4]a). 
As the index grows it is also shown that the contribution 
of short-period orbits become closer and closer to the one 
of the closed billiard. 

In Fig. [5] a) we present the spectrum of the TE res- 
onances for the square cavity of side a — 1 with ( ) 

symmetry along the diagonals. For such cavity the fit 
function similar to (j67l) is 

Nfitik) = —{kaf + ai{ka) + ao (70) 
167r 

and the best fit gives, see Fig. [5]b), 

ai = 0.0304, ao = -5.22 . (71) 



The theoretical prediction for this symmetry class is ob- 
tained ft-om ^ af" = 0.0297 which agrees well with the 
numerical calculations. 

Finally the same procedure was done for the dielec- 
tric stadium consisted of two half-circles of radius R con- 
nected by a rectangle with sides 2aR and 2R where a 
called the aspect ratio of the stadium. The calculations 
were restricted to the symmetry class such that the as- 
sociated function vanishes along both symmetry axis of 
the stadium, which is again called (— , — ) symmetry class. 
The resonance spectrum is presented in Fig. [5] a). 

The fit function is now 

2 

where the aspect ratio a has been taken to 1 in the nu- 
merical calculations. The best fit gives, see Fig. [HI b), 

ai = 0.150, ao = -5.24 , (73) 

which agrees well with the prediction for this symmetry 
class: af = 0.152 (cf. ([53])). 



VIII. SUMMARY 

Trace formulas are the main tool of the semiclassi- 
cal description of multi-dimensional quantum problems. 
For closed systems the trace formulas relate two objects: 
quantum density of discrete states and a sum over clas- 
sical periodic orbits 

11 

^d{E)+ Ape'^''(^)/^ + c.c. (74) 

period, orbits 

where Sp(E) is the classical action over a periodic orbit 
and d(E) is the mean density of eigen-energies, averaged 
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50 100 50 100 20 40 60 

kR kR kR 

a) b) c) 

FIG. 3. (Color online) Difference between the exact number of resonances and the fit (|67|l for n — 1.5 {&), n — 2 (b), and n = 3 
(c). In the latter case the additional resonances in Fig. [2]c) are not taken into account. The blue solid thick line indicates the 
difference averaged over a large interval. 



15r 




5.5 



1/R 




b) 



5.5 



1/R 



c) 



FIG. 4. Density of periodic orbit length (Fourier transform of the resonance density) for n — 1.5 (a), n = 2 (b), and n = 3 (c). 
The vertical lines stand for the length of the shortest periodic orbits of the circular cavity, from left to right: triangle, square, 
pentagon, hexagon, heptagon and octagon. 




FIG. 5. (Color online) a) Resonance spectrum for the dielectric square with n = 1.5 for the (— , — ) symmetry class, b) The 
difference between the total number of resonances and the best quadratic fit. The blue solid thick line indicates the difference 
averaged over a large interval. 
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Re(kR) kR 

a) b) 

FIG. 6. (Color online) a) Resonance spectrum for the dielectric stadium with n = 1.5 for the ( — , — ) symmetry class, b) The 
difference between the total number of resonances and the best quadratic fit. The blue solid thick line indicates the difference 
averaged over a large interval. 



over a small window around E. For 2d billiards with area 
A and perimeter C this averaged density of states is 



C 



+ o{E-'/^) 



(75) 



where r = 1 for the Neumann boundary conditions and 
r = — 1 for the Dirichlet ones. 

For open quantum models the true eigen-energy spec- 
trum is continuous and the main object of interest is the 
discrete spectrum of resonances defined as the poles of the 
5-matrix in the complex plane: E^ = e„ — ir„/2 with e„ 
and r„ real. The real part of the resonance energy, e„, 
gives the position of the resonance while its imaginary 
part, r„, determines the resonance width. 

The analogue of the trace formula for open systems has 
the form similar to (TM)) 



r„/2 



^Z^{E-e^Y+Tl/A 
« d{E) 



E 

period, orbits 



Ape'^p(iJ)/'« + c.c. 



(76) 



In Ref. [8| such type of formula has been obtained for 
a 2d dielectric cavity with transverse magnetic polariza- 
tion of the field. Here we derive the trace formula for 
a 2d dielectric cavity but with boundary conditions cor- 
responding to the transverse electric polarization of the 
electromagnetic field. As expected, the oscillating part 
of this trace formula is given by the usual periodic orbits 
weighted in the leading semiclassical order by the Fresnel 
coefficient corresponding to TE refiection on the cavity 
boundary (pS]) . (pSI) . 

Our main result is the expression for the average reso- 
nance density of a dielectric cavity with area A, perimeter 



C, and refraction index n 



d{E) = 4^ + rTE(n)-4= + o{E-''') (77) 



where 

rTE{n) = 1 + - / i?TE(i) 



l2 + i2 1 + t2 



± 



2n 



dt 



(78) 



and i?TE (i) is the Fresnel refiection coefficient for the TE 
polarization at imaginary momentum 



i?TE(t) 



i2 



(79) 



The plus-minus sign in front of the last term in (j78p is 
connected with the existence for the TE modes of an ad- 
ditional series of resonances related to Brewster's angle. 
As these resonances have large imaginary parts, they may 
be included or not in the counting function. For small 
values of n additional resonances are mixed with other 
resonances and their separation is artificial. In this case 
the plus sign has to be used. For large n the branch of 
additional resonances is well separated from the body of 
resonances and it is natural to ignore them. It corre- 
sponds to the minus sign in (1781) . 

The results of this paper together with Ref. [8| demon- 
strate that semiclassical trace formulas can be derived 
and applied for open dielectric cavities in a close similar- 
ity with closed billiards. Further investigations of trace 
formulas for other physical open systems is of consider- 
able interest. 
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Appendix A: Krein formula approach 

The purpose of this Appendix is to present another 
derivation of the number of resonances in a circular di- 
electric cavity based on the Krein spectral shift formula 
17| . The true eigen-energy spectrum for an open system 
is continuous and, consequently, the density of states for 
open quantum systems is infinite. Nevertheless, the dif- 
ference between the density of states with a cavity and 
the density of state without the cavity is finite and is 
given by the Krein formula 



d{E)-do{E) 



lndetS(£;) 



(Al) 



where S(i?) is the S'-matrix for the scattering on the cav- 
ity. 

This formula is general and can be used for any type 
of short-range potential. We apply it for a scattering 
on a circular dielectric cavity. It is easy to check that 
the S'-matrix for the the scattering on 2d circular cavity 
with TE boundary conditions ([5]) is diagonal in the polar 
coordinates and 



(A2) 



where Smix) is given by (O and Smix) differs from Sm(x) 
by changing Hm\x) to Hm\x): 



rl 



SM=x -J'^{nx)H^^Xx)-J-m{nx)H'i^'{x) . (A3) 
In 1 

From properties of the Bessel functions jisl it is straight- 
forward to show that 



'-{x) 



(A4) 



(n2 - 1) Jm{nx)Hl^\xW/x'' + nJ'^{nx)H^^' {x) 



J'^(nx)Hl^\x)ln - J^{nx)m7{x) 

Using the equality Jm{x) — {Hm\x) + Hm^)/2, this ex- 
pression can be rewritten in the form 



fkr N, _ n'^ - I Amix) + B„i{x)E„^{x) 



n2 Cm[x)\V - Rmix)Em{x)] 



(A5) 



where Em{x) and Rm{x) are defined in © and (flUl) re- 



spectively, and 

Am{x) 
Brn{x) 
Cm (x) 



m 

X 



H-(2)' rr(l)' 

11m iJ-m 
iiL 11m / \ -'-J-rn / \ 

_+n—^inx)-^ix) 

i^7n iJ-m 



(A6) 
(A7) 



H, 



(2)/ 



nH, 



■(nx) 



(1)' 



(A8) 



Expanding this expression into series of Em (x) gives 



■{x)^Qm{x)+PmY.R'm{x)E^:^{a 



where 
and 

Pm{x) 

with 



^ix) 



(n2 - l)Amix) 



n^Cm {x) 
4i(n2-l) TfS^ + (S^(-))' 



(A9) 



(AlO) 



'Kn?xHm\nx)Hm' {nx)Cm{x)Dm(x) 



r(2) 



(All) 



rr(l)' rr(l)' 

Dm{x)^^^{nx)-^{x) (A12) 

I^Elm Elm 

In the semiclassical limit x ^ oo the above formulae are 
simplified by using the asymptotic of the Hankcl function 

m 



H 



(1,2)' 



^0-.2) x^oo V x"^ 2{x'^ — m?) 
Consider first the smooth term (jAlOp . From the identity 
(^Vn2-t2 + VT~^) (y/n? - - VT~^) 

= + Vn2-<Vl-i2) (A14) 



ix) 



±i\ l- 



+ 0(x-2) . 

(A13) 



it is straightforward to check that 



1 - 



(A15) 



i?TE(^) 



where i?TE(i) is the Fresnel refiection coefficient for the 
TE polarization given by (IT5|) . 

The difference between the density of state with a cav- 
ity and the one without the cavity averaged over an en- 
ergy interval such that periodic orbit terms are small can 
be calculated from Qm{x) 



{d{E)}-do{E)^~^ IrnQmix) 



2-Kk 



(A16) 
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Changing the summation over m to the integration and 
turning the integration contour in the second term in 
(|A15P in the complex plane to avoid poles, m — > —it 
leads to 



(A17) 



{d{E)) - do{E) -- 


R - 


nnx j 

y f ^ 


27rfc . 


Rx 
A-Kk J 


' dtRTE 



—dm — 



1 Trdm 



(-4) 



1 



1^x2 +t'^ X^ + P 



Rescaling integration variables one gets 



{diE))^do{E) = ^in'-l) 



(A18) 



C 



dtRTsit) 



boundary conditions are well known (see e.g. |18| ) 



(2) 



(2)' 



(A20) 

The 'additional' S'-matrix for the TE scattering is thus 
formally 

— (m* + l) m* 

det(4™)(x))= n ^™^(^) n ^™^(^) 



m— — oo 

oo 



n ^^^^(^ 



(A21) 



m— m* + l 



where m* = [xn/\/ + 1], 

To find the total phase of this 'additional' S'-matrix 
one can proceed as follows. To the leading order in the 
semiclassical limit x ^ oo the Dirichlet and Neumann 
S- matrices (IA20p can be calculated from pT|) . It gives 



(A22) 



where A — ttR^ and C = 2ttR are the area and the 
perimeter of a circular cavity. 

This formula differs form the averaged total number of 
resonances (|4T|) and (|45l) . This is the consequence of the 
fact discussed in Q for the case of TM modes that the 
5-matrix for the scattering on a cavity has an additional 
phase (and additional zeros) connected with the outside 
scattering on the impenetrable cavity. 

The form of this 'additional' S'-matrix may be argued 
as follows. It is known that when a wave from outside 
the cavity scatters on a cavity it reflects with the reflec- 
tion coefficient which differs by the sign from the reflec- 
tion coefficient from inside the cavity (this is a conse- 
quence of current conservation). For the TE polarization 
the Fresnel reflection coefficient for a scattering from a 
medium with the refraction index 1 on another medium 
with the refraction index n is — i?TE where Rte is given 
by ((T5|) . In semiclassical region accessible in outside scat- 
tering, \t\ < 1, the reflection coefficient — i?TE(0 is real 
and the 'effective' reflection coefficient corresponding to 
the scattering on impenetrable cavity equals to the sign 
of -i?TE(t) (cf. dSSD) 



R. 



■TE 



-1, t* < |i| < 1 

1, < \t\ < t* 



(A19) 



where t* — nj \J"n? -|- 1. 

The reflection coefficient equals to —1 (resp. -1-1) cor- 
responding to the scattering with Dirichlet (resp. Neu- 
mann) conditions on the cavity boundary. For a circu- 
lar cavity the S-matrices with Dirichlet and Neumann 



where $,„ (2;) is given by . It means that Sm^ differs 
from Sm^ only by its sign, which is another manifestation 
of the opposite sign of the reflection coefficient (jA19l) . 
Therefore one can rewrite expression (|A2ip as follows 



,(TE) 



S(^) n (-1) « S(^)e±2iWV;i^ (A23) 



where S^-^^ = Vil=-oo ^^^^ is the fuU S-matrix for the 
scattering on a cavity with the Dirichlet boundary condi- 
tion. The ± sign in the exponent reflects the ambiguity 
of the phase, —1 = e^"^. 

The calculation of the mean density of states related 
with S'-^-'-matrix is straightforward (see e.g. [ist ) 



(A24) 



and finally from (|A23p and the Krein formula (|Aip one 
finds that the change of the density of states due to the 
'additional' S-niatrix (jA2ip is 



do(^) - do(^) = dD{E) - doiE) ± 

2n 



C 



2n 



87rfc 



A 



C 

8^ 



It 



1 



- 1 

(A25) 



The total density of resonances is thus the difference be- 
tween (|A18j) and (|A25p . In the end one gets Eqs. (|4T|) 
and (|45|) . The ambiguity in the phase of the 'additional' 
S-matrix corresponds to the the possibility to include 
resonances related with Brewster's angle in the Weyl for- 
mula or not which has been discussed in Section fVl 
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